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We study a spin S=M/2-Kondo system coupled to electrons in an arbitrary nonequilibrium situ- 
ation above Kondo temperature. Coupling to hot electrons leads to an increased inverse lifetime of 
pseudo particles, related to the Korringa width. This in turn is responsible for the increased inelastic 
relaxation rates of the electronic system. The rates are related to spin-spin correlation functions 
which are determined using a projection operator formalism. The results generalize recent findings 
for S=l/2-Kondo impurities which have been used to describe energy relaxation experiments in 
disordered mesoscopic wires. 
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Recently, experimental evidence was found that Kondo 
impurities might play an essential role for energy relax- 
ation in mesoscopic gold wires [l| displaying much higher 
energy relaxation rates than predicted by standard the- 
ory • Based on these findings several theoretical studies 
have led to a qualitative or even quantitative explanation 
of experimental data by accounting for electron-electron 
interaction mediated by magnetic impurities £1, Q, E| ■ 
Assuming Kondo impurities of unknown origin as rele- 
vant inelastic scattering centers also earlier experimental 
findings on copper wires could be explained 0, 0, @ • 
Moreover, assuming spin 1/2-impurities the detailed 
magnetic field dependence of energy relaxation experi- 
ments on copper wires could be fitted, strongly sug- 
gesting that Kondo impurities indeed play an essential 
role for energy relaxation at low temperatures. 

In a recent work Anthore et al. |lf)j | studied energy re- 
laxation in thin silver wires with Mn impurities and ex- 
plained their findings using both, direct electron-electron 
interaction 0] and the effect of spin 1/2-impurities 
Since Mn in silver is not a spin 1/2-impurity and further- 
more the spin of the impurities in copper is not known, 
a generalization of the theory in Ref. [8j is desirable. In 
addition the impurity densities c- lmp gained by fitting the 
energy relaxation data of the copper and silver samples 
typically exceed those obtained from measurements of 
the dephasing rate by more than an order of magnitude, 
see Refs. |3,ll0| an d articles cited therein. Impurity den- 
sities as high as those inferred from energy relaxation 
rates would lead to much higher dephasing rates than 
those found in experiments. 

Considering the theoretical work in Refs. 0, 0; the 
impurity density can be lowered by increasing the spin 5 
because only the product 5(5 + l)ci mp enters the prefac- 
tor of the rate. However, this result does not take into 
account the spin dependence of the renormalized coupling 
constant. The aim of the present work is a generalization 
of the findings in Ref. ;8<] to arbitrary spin thereby explor- 
ing the possibilities of lowering the impurity density by 
increasing the spin 5. 

Since this work is an extension of Ref. [8j we follow 
the argumentation therein and, as far as possible, use 



the same notation. In order to make the paper self- 
contained, some of the basic ideas and definitions are, 
however, repeated. Whereas the technical details change, 
the main physical arguments remain the same and we 
refer the reader to Ref. @ for further information. 

We describe the quasiparticles and the impurity spin 
by the free Hamiltonian 

H o = *^2 e k<r C la C kcr - E H S Z (1) 

ko- 

where c£ ff and C kcr create and annihilate an electron in 
a given orbital, k, and spin, a, state, eko- is the energy 
of this state. The second term in Eq. (JIJ describes a 
spin M/2-impurity with Zeeman splitting Eh = g^isB. 
The coupling between quasiparticles and impurity spin is 
described by the standard Kondo Hamiltonian 

H I = J S-s^Ct CT ,C kCT (2) 

kk'crcr' 

where Jo is the bare coupling and s a i a denotes the vector 
of Pauli matrices. Here, we assume the impurity density 
Cimp to be small enough that we need to treat coupling 
to a single impurity only. 

To determine the inelastic electron rates we consider 
the angularly averaged collision integral which in linear 
order in the density Ci mp reads 

I„{e) = ^|/ CT (e)E>(e) + [1 - /«r(e)]E<(e)j . (3) 

Here, E^^e) = S > / < (kcr, e) where e = is the elec- 
tron self-energy on shell, assumed to be independent of 
the angular momentum. f a (t) is the angularly averaged 
distribution function for electrons of energy e and spin 
a. For readability we suppressed the spatial dependence. 
Since the self-energy is proportional to the impurity den- 
sity, we already replaced the electron Green's functions 
by their unperturbed form and integrated over frequency 
to get the classical form of the collision integral. In con- 
trast to Ref. Q we do not use the spin averaged self- 
energy but generalize the results to spin dependent dis- 
tribution functions. 
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Our task is now to determine the electron self-energy 
which in turn leads to the electron scattering rates. 
Changing to the interaction picture and representing the 
spin degrees of freedom by pseudo-particles |T^ one can 
use perturbation theory on the Keldysh contour to gen- 
erate the graphs contributing to the electron self-energy. 
Since the topological structure of the graphs for a spin 
1/2-system and a spin M/2-system is the same, we can 
directly follow the reasoning in Ref. @- 

In lowest order the electron self-energy is given by a 
pseudo-Fermion bubble and an electron or hole line in be- 
tween. The pseudo-Fermion bubble can be represented 
as a spin-spin correlation function which in frequency 
space directly determines the rates. Higher order cor- 
rections are separable into terms adding an additional 
electron-hole pair and terms leading to higher order cor- 
rections for a single electron-hole pair. The combinations 
of the second type are usually referred to as singe par- 
ticle intermediate state corrections and can be absorbed 
by a renormalization of the coupling constants |13|. For 
arbitrary spin S, we find that the renormalized vertices 
J± , only depend on electronic occupation factors and 
therefore are given by relations very similar to those de- 
rived in Ref. y| . For a non-spin-flip processes we have 



J|(e)/J = { \l-(7TpJ ) 2 S(S+l)/4-pJog T (eTE H )\ 2 
+ {n P Jo?S{S+l)y l/2 (4) 
and for a spin-flip process 

J^/Jo = {\l-(7rpJo) 2 S(S+l)/4-pJ [g ± (e) 

+g T (e ± E H )}/2\ 2 + (tt P J ) 2 S(S + l)}~ 1/2 (5) 

The renormalization is determined by the auxiliary func- 
tion 



5±(e) = / de 

J-D 



,f±(e')-l/2 



iS 



(0) 



In equilibrium this leads to the usual logarithmic cor- 
rections, however, the above formulae are applicable for 
arbitrary nonequilibrium situations. The Kondo temper- 
ature in this approximation reads 



Tk = D exp \ — l — 
pJo 



1 - 



(n P J ) 2 S(S + 1) 



(7) 



and equals the bulk Kondo temperature. The phrase 
"above Kondo temperature" in this work means that the 
corrections determined by the auxiliary function © are 
still small compared to one. In this sense a system below 
the equilibrium Kondo temperature can be "above Kondo 
temperature" because of the nonequilibrium smearing of 
the distribution function. 

Well above Kondo temperature it is usually assumed 
that all vertices renormalize independently. Therefore, 



one can equivalently put these renormalized quantities 
in a new interaction Hamiltonian 

Irlr/ ^ 



+S Z [ji(e kT )Ct T C kT - J^(e ki )C k/i C ki 



(8) 



with energy and process dependent coupling constants. 
Using this Hamiltonian we have to restrict to elementary 
electron-hole pair excitations only. Other, more complex 
graphs of the one-particle intermediate state correction 
type, are already put into the renormalization of the cou- 
pling constants. The electron self-energy is now given by 
the pseudo-Fermion bubble coupled to arbitrarily many 
simple electron-hole pairs with an electron or hole line 
in between and can be written as 



S>(e) 



de'W a ^(e,e')[l-f a ,(e')} (9) 



for the larger self-energy where W a ,a< denotes the corre- 
sponding rates. The smaller self-energy £^(e) is given 
by changing the variables, (e, a) to (e', a') and / — > 1 — /. 
Rewriting the pseudo-Fermion bubble as spin-spin cor- 
relation function, the rates are given by 

CimpP 



W-,+ (e,e') = 
W+,-(e,e') = 
W+,+ (e,e') = 



m -J-(e)J + (6')C + (e- e ') (10) 

^J+(e)J-(e')C_(e-e') (11) 

^j;( e )4(e')C4e-e') (12) 

W_,_(e,e') = ^Pji( e )r(e')C z (e-e') (13) 

with 

C±(i) = (^(4)^(0)) , C z (t) = (S z (t)S z (0)) . (14) 

Using © the collision integral takes the standard form 
for spin dependent scattering 



E 



&'{/ CT (e)[l-/ CT ,(e')]^(e,e') 
-[l-/ CT ( e )]MeOW- CTV ( e ', e )}. (15) 



The energies, e = eko-, measure the kinetic energy and 
the Zeeman energy. Usually, when going over into a con- 
tinuum description the Zeeman splitting is put to a band 
bottom shift. 

As in Ref. [8j we use a projection operator formalism 
to determine the correlation functions in an arbitrary 
nonequilibrium situation. Using the projection operators 



P Z X 
P ± X 



S Z (XS Z )/(S Z S Z ) 
S ± (XS T )/{S ± ST) 



for C z and (16) 
for C± (17) 



one can derive a formally exact intcgro differential equa- 
tion m 

C a {t) = <5> a C a {t) - I du<t> a {t - u)C a {u) (18) 



3 



with the solution in terms of the Laplace transform 



C a (z) 



Cg(t = 0) 

Z - $ Q + 4> a (z) 



(19) 



Here a = z, ± and, as in the 5 = 1/2 case, $ z = 
(S Z S Z )/{S Z S Z ) = and $± = (S ± S^)/(S ± S^) = 
^fiEtf, which leads to the free propagation, where Eh in- 
cludes the Knight shift neglected throughout this work. 
The averages are to be calculated self-consistently to- 
gether with the steady state electronic distribution func- 
tions f a and the occupation probabilities V m for the im- 
purity spin being in state to. 

The memory kernel <fi a (t) for the C± correlation func- 
tion reads 



0±(t) = 



{S±S*) 



(20) 



Here, the index r in S^(t) indicates that the dynamics of 
the spin operator is reduced by the projection. It is de- 
termined by the expression S^(t) — exp[iL(l — P^iJS^ 
with the Liouville operator L acting as LX = [H, X] / h. 
The memory kernel for the C z correlation function is 
given by Eq. 1)20(1 with the replacements ±, T - * z - We 
are interested in the regime well above Kondo temper- 
ature and expand the kernel up to second order in the 
renormalized coupling J. Since the dynamics of the ex- 
panded kernel function is oscillatory, the Fourier trans- 
formed correlation function has always the simple form 



Ca(Lu) 



2C a {t = O)Re0„(w) 



[u - i$ a + Im^H] 2 + [Re^ a (w)] 5 



(21) 



with a = z, ±. Further, we define Re0 a (cj) = 
Re{0 a (— iui + 5)} and the imaginary part Im<^ a (w) fol- 
lows from the Kramers-Kronig relation. When calcu- 
lating the electronic distributions f a or spin occupation 
probabilities V m , the imaginary parts Im0 a (o;) in the de- 
nominators lead to higher order corrections in J and are 
neglected. 

The damping rates (which were named v a (oj) with a = 
z,± in Ref. @) read 



(5±5^) 



(S Z S Z 



Q ± {ujtE h ) 



for the C z correlation function and 



Re0±(w) = - 



iqz cz\ 



(22) 



(23) 



~(5±5T) 

for the C± correlation functions. The auxiliary functions 

CH = E / ^ 2 J|(e)J|(e + ^)/ ± (e)[l-/ ± (e + c)] 
± J 

(24) 



and 

C±(w) = ydep 2 jT( e )J±( e -|- w )/ T ( e )[l-/ ± ( e + w )] (25) 

describe coupling to electron-hole pairs. In equilibrium 
the damping leads directly to the Korringa width propor- 
tional to the temperature whereas in noncquilibrium this 
rate scales with a measure of the nonequilibrium situa- 
tion, namely eU, leading to an increased inverse lifetime 
independent of the measurement temperature. 

The equal time correlation functions read for 5 = M/2 

M/2 

C z (t = 0) = (S z S z )= £ V m m 2 (26) 

m=-M/2 



and 



M/2 



C ± (t = 0) = {S ± S*)= ]T V m [S(S + l)-m(mTl)}. 

m=-M/2 

(27) 

Independent of the distribution V m the spin-spin correla- 
tion function C(t) = (S(i)-S) - [C+{t) + C-{t)]/2+C z (t) 
fulfills the sum rule C(t = 0) = f(du)/2n) C(u>) = 
5(5 + 1). 

To determine the master equation for the P m 's we use 
Eq. 127|) and write the spin-flip correlation function as 

OtM = E P m [S(S + 1) - m(m T 1)]<5±M . (28) 

m 

The rate for the transition from state to to to ± 1 then 
reads 

r m ^ m±1 = [5(5+l)-m(m±l)]r± (29) 

with 



r ± = ^ / dwC=F(-w)GFM. 



(30) 



All other rates vanish. Note, that the definition of T± in 
this work is different from that employed in Ref. . The 
rate equations for the occupation probabilities 

r m = —T m ^ m j r iP m — r m _ >m _i'P m 

+F m +l^m'P m+1 + T rn -i^ m 'P m _ 1 (31) 

with the normalization condition ^2 m V m = 1 form a 
closed set of equations with the steady state solution 



V. = 



,M /2+m r M/2-m 
+ 1 - 



EJw r-\M — n-pn 
n=0 i + - 



(32) 



The probabilities obey the obvious balance relation 
V m /V m+ i = r_/r + which leads to the thermal distri- 
bution in equilibrium. 

At vanishing magnetic field, B = 0, the probabil- 
ities are all equal, V m = l/C^ + 1)) an d the equal 



4 



time correlation functions read C z (t = 0) = 5(5 + l)/3 
and C±(t = 0) = 2C z (t = 0). If in addition the 
distribution functions are spin independent, the renor- 
malized coupling constants become process independent 
J± = = J(e), and the auxiliary functions read Q z — 
2(± = 2((uj). Inserting this in the correlation functions, 
we find C(u) = [C + {lo) + C-{lo)]/2 + C z {lo) = 3C+(w)/2. 
In equilibrium and at low temperatures the width shrinks 
to zero and leads to C(ui) —* 2nS(S + l)S(w). 

The inelastic relaxation rate l/T; no i at B = is the 
spin-flip rate l/r sf reduced by the quasi-elastic rate, and 
in general we have l/T- mc \ < l/r s f. Quite generally, due 
to a sum rule for the spin-spin correlation function, the 
spin-flip rate obeys 

Tet 2 ^ J 2h p 

a, a 

In order to discuss the possibility of reducing the im- 
purity density by increasing the spin 5 at constant in- 
elastic electronic rate we may as well consider the spin- 
flip rate. As already explained in the introduction, Eq. 

suggests a decrease of the impurity density with in- 
creasing spin 5. This is true only if the renormalization 
of the coupling constants is independent of 5 meaning 
at temperatures much higher than the Kondo tempera- 
ture. To explain the experiments, however, pJ has to be 
around 1/3 to be almost voltage independent. Otherwise 
the renormalization wouldn't allow for the experimen- 
tally observed scaling property of the distribution func- 
tion f(e,eU) — f(e/eU), see Ref. 0,0. In this regime, 
however, the renormalization depends on the spin 5 and 
scales for large spin like pJ ~ l/^n 2 S(S + 1) leading to 
a spin independent rate l/r s f. Actually, the renormal- 
ized coupling constant equals the spin-flip t-matrix p| 
which obeys a unitarity condition. It reaches a maxi- 
mum, p J = 1/ y / Tr 2 S(S + 1), at the Kondo temperature 
where the rate again would become independent of spin 



for all 5. Although our theory is no longer valid in this 
regime, the outcome is quite physical since electrons al- 
ways transfer the same spin when scattering from one 
impurity independent of 5. This shows that using our 
theory an increase of the spin does not lower the impu- 
rity concentrations needed to describe the experiments. 
Even a more involved theory valid below Kondo temper- 
ature is not likely to help much since the scattering rate 
cannot exceed the limit discussed above. 

To discuss the magnetic field dependence of the rates 
W we consider two limiting cases. For low magnetic fields 
where the Zeeman splitting Eh is much smaller than the 
temperature or the applied voltage the occupation prob- 
abilities are all of the same order. Also the lifetimes do 
not change much and the behavior is dominated by the 
shift in the spin-flip correlation functions, u> — > U) ± Eh- 
Therefore, there is no dependence on the spin 5 for small 
magnetic fields. For higher magnetic fields of the order 
of temperature or applied voltage, higher spin states are 
rapidly depopulated so that only two spin states like in 
the 5 = 1/2 case lead to the dominant contribution. For 
higher 5 this is of course just a fraction and therefore in 
this regime the rates are even smaller than in the 5=1/2 
case. 

In this work we have studied electron relaxation rates 
caused by magnetic impurities of arbitrary spin gener- 
alizing recent results for 5 = 1/2. It is found that an 
increase of the spin 5 does not change the qualitative 
outcome and the rate at vanishing magnetic field is even 
unaffected by the spin for large 5. Therefore, assuming 
magnetic impurities with higher spin 5 does not resolve 
the disagreement between Kondo impurity densities de- 
termined by energy relaxation experiments and weak lo- 
calization experiments. The authors would like to thank 
B. L. Altshuler, A. Anthore, Y. M. Galperin, F. Pierre, 
and H. Pothier for valuable discussions. Financial sup- 
port was provided by the Deutsche Forschungsgcmcin- 
schaft (DFG). 
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